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ABSTRACT 

We present new predictions for temperature (on small angular scales) and polarization (on large 
angular scales) CMB anisotropics induced during the epoch of reionization (EoR). Using a novel 
method calibrated from radiation-hydrodynamic simulations we model the EoR in large volumes (L 
> 2 Gpc/h). We find that the EoR contribution to the kinetic Sunyaev- Zel'dovich power spectrum 
(patchy kSZ) ranges between ^0.6 - 2.8 ^K 2 at I = 3000, for the explored parameter space. For each 
model, the patchy kSZ power spectrum is calculated from 3 large 15° x 15° maps for better numerical 
convergence. Decreasing the size of these maps biases the overall patchy kSZ power to higher values. 
We find that the amplitude of the patchy kSZ power spectrum at I = 3000 follows simple scalings of 
-Dfciooo ^ ^ an d -Dfclooo ^ ^z' 47 for the mean redshift (z) and duration (A z ) of reionization. Using 
the constraints on z from the WMAP 7- year results and the lower limit on A z from EDGES wc find a 
lower limit of ~ 0.4 fjK 2 at I = 3000. Planck will constrain the mean redshift and the Thomson optical 
depth from the low-£ polarization power spectrum. Future measurements of the high-^ CMB power 
spectrum from the South Pole Telescope (SPT) and the Atacama Cosmology Telescope (ACT) should 
detect the patchy kSZ signal if the cross correlation between the cosmic infrared background and the 
thermal Sunyaev Zel'dovich effect is constrained. We show that the combination of temperature and 
polarization measurements constrains both z and A z . The patchy kSZ maps, power spectra templates 
and the polarization power spectra will be publicly available. 

Subject headings: Cosmic Microwave Background — Cosmology: Theory — Galaxies: Clusters: Gen- 
eral — Large-Scale Structure of Universe — Methods: Numerical 



1. INTRODUCTION 

Free electrons in the intcrgalactic medium (IGM) scat- 
ter photons from cosmic microwave background (CMB) 
creating additional secondary anisotropics that distort 
the primordial anisotropics. These free electrons are ini- 
tially ionized from the neutral IGM by the first stars and 
galaxies during the epoch of reionization (EoR). Thus, 
information on EoR is imprinted on the CMB in both 
temperature and polarization. The temperature fluctua- 
tions are affected by the kinetic Sunyaev Zel'dovich (kSZ) 
effect, which is doppler shifting of CMB photons from the 
bulk motio ns of free electrons with res pect to the CMB 
rest frame ([Sunyaev fc Zeldovichl[L980| ).The polarization 
signal on large angular scales is affected by the rotation 
sourced by free e lectrons from the beginning of EoR to 
the present ('e.g. iBond fc Efstathioul 119871: IHu fc White 
1997b iZaldarriaeal 119971: iLiu et all 120011: iKomatsu eTaL 
20111) . which induces a curl free polarization signal (E- 
mode). These are the two leading order effects, while 
there are other small er order effects, such a s fluctuations 
in the optical depth (|Natarajan et al.ll2012[ ) , that are not 
discussed. 

Alrea dy measurements from the opacity of the Lya 
fores t (lEan et all l2006bl ) the redshifted 21 cm sig- 
nal (| Bowman fc Rogers! I2012T) from the experiment 
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EDGES 4 , and the l arge-scale polarization of the CMB 
(|Larson et al.l 1201 lh infer that reionization was ex- 
tended. Further model dependent constraints on EoR 
come from: measurements of quasar prox imity zones 
(e.g. iWvithe et ahl [20051 iFan et all [20061 . a null re- 
sult for intergalactic damping wing abs orption in a 
z = 6 .3 gamma-ray burst s pectrum (e.g. iTotani et all 
120061: iMcQuinn etall 120081) . detections of damping 
wing absorption in the IGM from quasar spectra (e.g . 
Mesinger &: Haimanl [2004: Mesinger fc Furlanettcl 12001 



Bolton et al.l 120111 ). and Lya emit ter number densities 



and c lustering measuremen ts (e.g. iMalhotra fc Rhoadsl 
l200i iHldnTan fc Cenl [20051) . Recently, the South Pole 
Telescope (SPT) placed a model dependent upper limit 
on the duration of reionization (jReichardt et all 120121 : 
IZahn et"aTI 120121 ) from their multifrequency measure- 
ments of the high-£ power spectrum of CMB secondary 
anisotropics. New CMB measurements of temperature 
and polarization anisotropics from the Planck satellite, 
the POLARBEAR experiment, the Atacama Cosmol- 
ogy Telescope (ACT), ACT-pol (ACT with polarization), 
SPT, SPT-pol (SPT w ith polarization), and CMBpol 
(jZaldarriaga et alll2008l ) have the potential to constrain 
the EoR from CMB measurements alone. 

The amplitude of the EE power spectrum at £ < 20 
essentially measures the optical depth to reionization, r, 
with the most recent constraints being t = 0.087 ±0.015 
and a mean reionization-redshift is 10.5 ± 1. 2 (68% CL). 
This amplitude is often predicted by CAMB ([Lewis et all 
120001 ) . or codes like it, using a parametric hyperbolic 
tangent function for the ionization history. There are 
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modifications to CAMB (jMortonson fc Hul 120081 ) . which 
allow any ionization history as an input. 

The fract ional contributi o ns to the kSZ from the 
EoR (e.g. iGruzinov fc Hul 119981: iKnox et al.1 11998 



Valaeeas et alJ 2001: Sa ntos et al . 2003; Zahn et al .ll2005l 
McQuinn et all 120051: llliev et all 120071: IMesinger et al.l 
20111: iZahn et al.l 120121 : iMesinger et"afl 120121 hereafter 



we refer to this contribution as the patchy kSZ) are 
the largest on small angular scales compared to the 
primary and other secondary CMB anisotropies. This 
patchy kSZ power is in addition to the kSZ power 
that c omes from lower reds h ift (|Ostriker fc Vishniad 
19861 Uaffe fc Kamionkowskl Il998t IMa fc Frvl 120021 



Zhang et al.l 120041: iShaw et al.l I2012L here after we refer 
to this contribution as the homogeneous kSZ). Many of 
the previous models for the patchy kSZ signal were calcu- 
lated in small volumes (< 1 Gpc/h) and do not capture 
the large scale features of the patchy kSZ maps, which is 
required to accurately calculate the power spectrum. 

This is the third paper (Paper III) in a series that ex- 
plores EoR observables produced via our semi-analytical 
models of reionization that are statistically informed 
by simulations with radiative transfer a nd hydrodynam- 
ics. W e introduce our model in Paper I (|Battaglia et al.l 
l2012bf ) . we look at the impact of a patchy optical 
depth on CMB observables in Paper II (|Nataraian et al.l 
120 12f). and we explor e the 21cm signal in Paper IV 
(|La Plante eta[MoT% . 

We present in this paper predictions for CMB observ- 
ables. These predictions are made in large volumes (L 
= 2 Gpc/h) and the importance of going to such large 
volumes is demonstrated throughout this work. The EE 
polarization power spectra, the kSZ power spectra, and 
the maps from this paper will be made publicly avail- 
able. In Section[2j we summarize our fast semi-analytical 
model and the simulations it is based on. In Section 02 
we present results for the EE power spectrum and the 
kSZ power spectrum. We discuss prospects for future 
measurements and conclude in section [5j We adopt 
the concordance cosmological para meters that are con - 
sistent with WMAP 7-year results (jLarson et al.l 1201 If ): 
n m = 0.27, n A = 0.73, n h = 0.045, h = 0.7, n s = 0.96, 
and (j g = 0.80. 

2. PARAMETRIC MODEL FOR REIONIZATION 

In Paper I we developed a semi-analytic model 
for reionization based upon results from RadHydro 
simulations (for mor e details see iTrac et al.l 120081: 
iBattaglia et al.|[2012bf ). In the simulations we construct 
a reionization- redshift field, zre(x), that tracks the red- 
shift at which each the gas cell becomes 90% ionized. We 
define the following fluctuations fields for density 



^m(x) 



P(x) 



p 



and the reionization-rcdshift 

[1 + z RE (x)] - [1 + z] 



4(x) 



1 + z 



(1) 



(2) 



where p is the mean matter density and z is the mean 
value for the zre(x) field, which is approximately equal to 
the redshift of 50% ionization. The fluctuations in both 
the <5 m and S z fields are highly correlated on scales > 1 



Mpc//i. We calculate a simple scale-dependent linear 
bias that relates these two fields and we represent this 
bias with the simple parametric form, 



b zm (k) ■■ 



(S z (k)S z (k)) 



{S m (k)S m (k)) 

(l+fc/fco)"' 



1/2 



(3) 



that contains 3 parameters b a , k a , and a. The value for b 
that we use is de t ermin ed from analytical arguments in 
iBarkana fc Loebl (|2004l ). The fiducial parameter values 
for k — 0.185 Mpc/h and a — 0.564 are found by fit- 
ting the bias calculated from the simulations. We explore 
the parameter space of our model by varying k a and a. 
The effects these parameters have on EoR are the follow- 
ing: increasing k Q lengthens the duratio n of reionization 
while i ncreasing a shortens reioinization (Battagli a et al.1 
I2012bf ). Physically, a shorter reionization process tends 
to have larger ionization bubble sizes that percolate more 
quickly, which in turn correspond to more luminous ion- 
izing sources. 

We generate the over-density fields, S m using a particle- 
particle-particle-mesh (P 3 M) N-body code that evolves 
2048 3 dark matter particles in a 2 Gpc/h box down to 
z = 5.5. This over-density field is then convolved with 
a filter consisting of three elements: (1) a cubical top 
hat filter, which deconvolves the smoothing used to 

construct S m from the simulation, (2) a Fourier transform 
of a real space top hat filter 9 (k), which smoothes <5 m to 
a resolution of 1 Mpc/h, and (3) the bias function from 
Equation [3] The assembled filter takes this form 



W z {k) = 



b zm (k)Q(k) 



(4) 



and we apply this filter at z. The newly constructed 
S z field is Fourier transformed back to real space and 
converted to the zre(x) field by Eq. [5] with the same z 
as the density field. Here the value of z essentially sets 
the midpoint of reionization. We now have a complete 
ionization history for the density field used, which is then 
used to make ionization fields and kSZ maps. We define 
the duration of reionization as 



z(x c = 25%) -z(x e = 75%), 



(5) 



where x a is the ionization history. This definition ex- 
cludes the early and late times of reionization, since the 
small scale physical processes at these times are difficult 
to capture, hence is insensitive to the epochs when both 
our simulations and semi-analytic method are most un- 
certain. For a detailed parameter study of A z (fc Q , a) see 
IBattaglia efall (|2012bD . 

3. RESULTS 

We present our model predictions for the integrated 
optical depth, r, the low-^ EE mode polarization power 
spectrum, and the contribution to the patchy kSZ power 
spectrum. These predictions are compared to current 
and projected constraints, as well as previous work. 
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Fig. 1. — Left: The ionization fraction as a function of redshift, x e (z), for five models, our fiducial model at at z = 8,10, and 12 (green, 
red, and blue, respectively) and two extreme models of brief (orange dashed) and long (light green dashed) duration rcionization at a fixed 
z = 10. The x e (z) for the fiducial models have similar shapes and they are just shifted according to z. Right: The correspo nding optical 
depth , r, for the same models as x e . The values of r are compared against the WMAP 7-year constraints (light grey box Larson ct al. 
120111) on t and their reionizatoin-redshift. Like the polarization power spectrum, a constraint on r will differentiate between models with 
different z, but cannot help differentiate between our models with large or small A z . 
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Fig. 2. — Top: The \ow-£ EE polarization power spectrum 
for several models of reionization compared to the WMAP 7 - 
year band power constraint (light grey box [Larson ct af] 120111 ) 
and the projected error bars from Planck (dark grey bands 
IThe Planck Collaboration! 12009 ) for the 143 GHz channel calcu- 
lated for our fiducial model. We show five models, our fiducial 
model at at z = 8,10, and 12 (green, red, and blue, respectively) 
and two extreme models of brief (orange dashed) and long (light 
green dashed) duration reionization at a fixed z = 10. Bottom: 
The percent difference between our fiducial and the two extreme 
models. The low I EE polarization power spectrum will tightly con- 
strain z. However, it is not possible to discern between our models 
with large or small A z for same z using only this measurement. 



.1. Optical Depth and EE Polarization Power 
Spectrum 



We calculate the lar ge scale EE pola rization power 
spectrum us ing CAMB (|Lewis et al" I2000T) with the mod- 
ifications bv lMortonson fc Hul (|2008 ). which accepts gen- 
eral reionization histories, x e {z). Although novel, we 
did not include their PCA implementation. Figure [T] 
shows the results for x c (z) from our fiducial models at 
z = 8, 10, 12 and two extreme models for long and short 
duration reionization scenarios at z = 10. We show the 
corresponding EE polarization power spectra in Figure 
[2] and comp are them to the W MAP 7-year band power 
constraints ([Larson et al.ll20iH) . We find the fiducial pa- 
rameters with z = 10 agree with the WMAP constraint, 
and the 1-a confidence interval is bracketed by choices 
of z = 8 and z = 12. Our two extreme models of short 
and long duration reionization (here after small and large 
A z ) at fixed z differ from the fiducial model by a max- 
imum 20% (cf. Fig [2). Illustrated in Fig. |2] is that 
the current WMAP 7-year data is unable to differen- 
tiate between these extreme models for reionization and 
given the projected error bars on the EE power spectrum 
(|The Planck 'Collaborationll2006l for the 143 GHz chan- 
nel and our fiducial model) neither will the upcoming 
observations from Planck. Thus, the low £ EE polariza- 
tion power spectrum constraints are not sensitive to the 
duration of reion ization. This conc lusion is the same as 
previous work bv lZahn et al.l (|2012f) . however, they come 
to this conclusion via a different semi-analytic model. 
Similarly, WMAP showed that only using primary CMB 
constraints there is a degeneracy betwe en the duration 
and the mean redshift of rcionization dKomatsu et al.l 
[20ll . There is no benefit to include the low £ TE cross 
spectrum, since the projected e rror bars from Planck 
(|The Planck Collaborationll2006f ) are too large to differ- 
entiate between small and large A z models. 

Any constraint on r is primarily driven by the mea- 
surement of the low £ EE polarization power spectrum. 
To first order the amplitude of the low £ EE power spec- 
trum goes like r 2 . In Figure Q] we compare the WMAP 
7-year results to our results for the integrated r which is 
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Fig. 3. — Light cone projected maps of the kSZ signal from patchy reionization for models with z = 10 and our longest A z (left), our 
fiducial A z (center), and our shortest A z (right). The overall large-scale structure is similar but the small-scale structure decreases as A z 
decreases. 
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Fig. 4. — The patchy kSZ power spec trum for various r eionization 
models compared to SPT constraints (Zahn ct al. 2012). We show 
five models, our fiducial model at z = 8,10, and 12 (green, red, and 
blue, respectively) and two extreme models of brief (orange dashed) 
and long (light green dashed) duration reionization at a fixed z = 
10. The SPT constraints are illustrated by the grey arrows, with 
the darker grey arrow representing the constraint ignoring the tSZ- 
CIB cross-correlation and light including this correlation. All our 
models fall below the constraint that allows for tSZ-CIB correlation 
with an I dependent shape and most of the models are below the 
even tighter constraint which ignores this correlation. 



a consistency check on the EE power spectrum results, 
since there is no extra information. Comparing Figures [T] 
and [2] we find that the models with larger r arc consis- 
tent with having larger EE power at low-£ Like the EE 
power spectrum, measuring r places a constraint on z 
but is unable to discern between models with small and 
large A z . 

Here the assumption is that r is uniform in all direc- 
tions, however, reionization is naturally inhomogeneous. 
The optical depth as a function of the distance from an 
observer to the CMB is given by the equation, 



T(l*,n) = (JT I n e (n,l)dl, 
lo 



(6) 



where ctt is the Thomson cross section, n is the direc- 
tion normal unit vector, dl — cdt is the proper distance 
along the line of sight, U is the distance to the surface of 



last scattering, and n c (h,l) is the ionized election num- 
ber density at position I in the direction n. Previously, 
there was an upper limit constraint put on the RMS fluc- 
tuation in r of, at most, a few p ercent of the mean valu e 
(t) using published SPT data (jMortonson fc Hull20ld ). 
Given this upper limit the patchy r contribution to the 
CMB power s pectrum at high-l is n egligible. In a com- 
panion paper (jNataraian et al.ll2012f ). we show that four 
point statistics of the CMB, in principle, can constrain 
the RMS fluctuation in r and if measured one can dif- 
ferentiate between models with small and large A z and 
possibly break this degeneracy between A z and z . 

3.2. Patchy Kinetic Sunyaev ZeVdovich Power 
Spectrum 

The kSZ signal from patchy reionization is sensitive to 
the details of reionization such as the mean redshift of 
reionization and its duration (e.g. IMcQuinn et al.l I2005t 
IZahn et al.ll200l IMesinger et alj|2012tlZahn et al.ll2012D . 
which, physically, is primarily due to the relatively steep 
dependence on redshift of the signal strength and the de- 
pendence of ionization bubble size on duration. In this 
section, we investigate how the patchy kSZ power spec- 
trum depends on our model parameters and we compare 
our results against previous work and observational con- 
straints. We provide a simple scaling relation for the 
patchy kSZ power at I = 3000 as a function of z and 
A z , which makes model fitting of observational spectra 
trivial. 

Wc construct patchy kSZ maps by raytracing through 
the past light cone (5.5 < z < 20). The temperature 
distortion along each line of sight is given by 



AT 



■(n) 



or 

c 



r(i,: 



'n e (l, n)n • vdl, 



(7) 



where v is the peculiar velocity, r(n) is from Eq. [SJ 
and l is the proper distance at z = 5.5. We make flat 
sky maps that are approximately 15° x 15°, where the 
angular size is determined by the N-body simulation box 
length of 2 Gpc/h over the comoving distance out to 
z = 20. Since the box length (L = 2 Gpc/h) of the N- 
body simulations is approximately equal to the comoving 
distance between z = 6 and z = 20, wc cycle through 
the projection direction coordinates approximately once 
when making the maps. For each choice of parameters 
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Fig. 5. — Parameter study of the patchy kSZ power spectrum at 
I = 3000 for varying k a and a with a fixed z = 10. Th e black line 
indicates the upper limit from SPT (Zahn ct al. 2012) with £ = 
and all models above that are consistent with this constraint. The 
bottom right corner has the largest A z and the top left corner has 
the shortest A z . 



we make three maps along three independent axes. In 
Fig. [3] we show the patchy kSZ Compton-y maps for our 
fiducial model and the two extreme models for the same 
projection direction. From these maps it is obvious how 
the duration of ionization affects the patchy kSZ. The 
models of reionization with large A z have more small 
scale structure than models with small A z . 

Using the flat sky approximation we calculate the 
power spectrum from the patchy kSZ maps and average 
over each projection axis for a given parameterization. 
Figure U shows the patchy kSZ power spectrum for var- 
ious models of reionization from the 15° x 15° patches. 
We find that increasing z with a fixed A z increases the 
overall amplitude of the patchy kSZ power spectrum, but 
has little effect on the shape since by z docs not affect 
the shape of x c {z) or A z . Altering A z at fixed z dra- 
matically changes both the amplitude and shape of the 
patchy kSZ power spectrum, sin ce A 7 affects the correla - 
tion between ionization regions (Battagl ia et al.l[2012b[ ). 
Models with smaller A z have large ionization regions, 
and more power at smaller I compared to larger A z mod- 
els where the power peaks at larger I due to the smaller 
coherent ionized regions. These results are qualita tively 
similar to [Mesinger et all (|2012D ; lZahn et al.l p012h . The 
strong shape and amplitude dependencies of the patchy 
kSZ power spectrum on A z illustrates that if there are 
constraints on z from the EE power spectrum and r, then 
a measurement of the patchy kSZ power will break the 
degeneracy betwee n z and A z , which was already shown 
using SPT results (jZahn et al.ll2012T) . 

In Figure |4j our patchy kSZ powe r spectra are com - 
pared against the SPT upper limits (jZahn et al.ll2012l ). 
which accounted fo r the homogenous kSZ contribution 
using a model from lShaw et al.l (|2012l ). Here the tight- 
est constraint from SPT does not account for any cor- 
relation, £, between thermal SZ and the cosmic infrared 
background (CIB), which they measured at I = 3000 to 
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Fig. 6. — The percent difference between the average patchy kSZ 
power spectra from small cut out maps compared against our 15° X 
15° maps using the fiducial parameters with a z = 10. The average 
power from the 7.5° X 7.5° cut out maps is blue line and the average 
power from the 3.75° X 3.75° cut out maps is the red line. All these 
small maps have more power than our 15° X 15° maps. 



be &=30oo = -0.18 ± 0.2 (jReichardt et al.ll2012T) . The 
other two constraints account for a non-zero £, where £ 
is given a uniform prior from -1 to 1 and either includes 
i dependent shape constraint or does not. We find that 
our fiducial model is consistent with all SPT upper limits 
regardless of tSZ-CIB correlation treatment. Not all our 
models are consistent with SPT upper limits, for exam- 
ple models where we increase z or A z . Figure [5] shows 
our parameter space study of a and k for the patchy 
kSZ power at I = 3000 with a fixed z = 10. The ampli- 
tudes of the patchy kSZ power in Fig. Orange from 0.87 
- 2.42 /iK 2 and some models with large A z (i.e. low a 
and large A: Q ) do not fall below the SPT constraints that 
exclude the tSZ-CIB correlation (cf. Fig. S]). 

The kSZ power spectra that our model produces are 
lower than previous pr edictions from semi-analytic mod- 
els and simulatio ns (e.glMcQuinn et al .11200511 Zahn et al 
20051 llliev et al.ll2007t iZahn et al.ll2012t iMesinger et al 



20121 ). Some of these earlier works used much smaller 
volumes to calculate their patchy kSZ signa l 



iMcQuinn eFail [20051: IZahn et al.|[2005l: llliev et al" 

The later wor ks either used mo dels with restricted pa- 
rameter space (|Zahn et al.l l2~0121 to fit the data, or used 
the Limber approximatio n to calculate the pa tchy kSZ 
from 3D ionization fields (|Mesinger et al.ll"2012l ). Testing 
these assumptions is left for future work. 

We calculated the patchy kSZ using the flat sky ap- 
proximation on 15° x 15° maps. There are large struc- 
tures in these maps that would not be captured by small 
scale maps (cf. Fig. |3]). We checked how calculating the 
patchy kSZ power spectrum on smaller maps would bias 
the patchy kSZ signal by dividing our maps into four 
7.5° x 7.5° and sixteen 3.75° x 3.75° maps (totaling 12 
and 48 maps, respectively) and calculating the spectra 
for each new map and averaging them together. In Fig- 
ure [6l we show that compared to our original kSZ power 
spectrum the average power spectrum of the small maps 
had more power. At i = 3000 the fractional differences 
are ~ 2% and - 10% for the 7.5° x 7.5° and 3.75° x 3.75° 
maps, respectively. In the flat sky approximation small 
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Fig. 7. — The scaling relations between the patchy kSZ amplitude at I = 3000 and A z and z, here the red symbols are the values 
calculated from the maps, the blue lines are the best fit, and green symbols are the values from the combined scaling law (cf. Eq. llOI I. Left: 
The patchy kSZ amplitude at I = 3000 as a function of z scales like f^£|ooo ^ ^he squares, crosses, and triangles represent A z = 0.2, 
1.05, and 2.10, respectively. Right: The patchy kSZ amplitude at £ = 3000 as a function of A z scales like D^^ Q0Q oc Aj ,s , here z = 10. 
This combined scaling relation can be used to constrain EoR in secondary parameter fitting of the high-^ 1 CMB measurements. 



area maps are susceptible to erroneously produce more 
power on all scales since there are large scale features 
that these small maps do not capture. 

We found that the dependence of the amplitude of the 
patchy kSZ power spectrum at £ = 3000 on z and A z can 
be represented by simple scaling laws. Figure[7]illustrates 
these dependencies of £>^| 000 = C^ ma £(£ + l)/(2n), 
which is linear for z and a power law for A z . When fit- 
ting for the scaling laws we use a nonlinear least-squares 
method where each value is weighted by the inverse of the 
variance of the three different projections and treat the 
dependence of D*^ 000 on z and A z as separable func- 
tions. The scaling laws are constrained to be: 



for a fixed A z = 1.05 and 









-0.14 







A z 
1.05 



0.47 



(8) 



(9) 



for a fixed z = 10. The combination of Eq. [8] and [9] 



r> KS 

u e— 



kSZ 



£=3000 



2.02^K J 



11 



-0.12 



A z 
1.05 



0.47 



(10) 

The predicted D^E| 000 from Eq. [TO] compares well to 
the results from our model (cf. Fig. [7]). The slight devia- 
tions are seen in the variation of the map power spectrum 
values about the fit (right panel Fig. [7]) and these devia- 
tions are found at the extreme ends of parameter space. 
Using Eq. [TO] we find a lower limit of £>fc| 000 > 0.4/uK 2 
by taking the 2<r lower confidence interval on z = 8.1 
from WAMP7 and the lower limit on A z > 0.07 from 
EDGES. Here we have converted the EDGES definition 
of A z , which assumes a functional form of hyperbolic tan- 
gent for x e (z) to our definition. This scaling law provides 



a simple way to place model dependent constraints on z 
or A z by including it when fitting high-£ CMB power 
spectra measurements into the secondary models used. 
However this requires additional measurements, for ex- 
ample of the EE power spectrum, to break the degen- 
eracies between A z and z that occurs when just using 
patchy kSZ measurements. 

4. FUTURE CONSTRAINTS 

Constraints on z and A z will tighten as the preci- 
sion increases on measurements of the low-(? EE polar- 
ization and high-£ temperature power spectra. We fore- 
cast how well these future precision measurements of r 
and -D£ = 3ooo will constrain z and A z by constructing a 
likelihood surface from a x 2 grid of z and A z around our 
fiducial model. The x 2 grid is calculated following 



x 2 = 



T - Tfid 


2 


^£=3000 — -C>£=3000,fid 


2 




+ 


^-D^ = 3000 





(11) 



here r^d and £^=3000, fid are the values for t and -Df = 3ooo 
from the fiducial model, cr T is the forecasted error bar for 
Planck or CMBpol on r, and &D e=3000 is the hypothetical 
error bar for ACT-pol and SPT-pol on Di =30 oo- Using 
our fiducial value for r we esti mate that Planck will mea- 
sure t ± 0.004 (~ 5% error; iThe Planck Collaboration] 
[2006] ) and C MBpol will measure r ± 0.002 a factor 
of 2 better (|Zaldarriaga et all 120081) . There is still 
no detection of the patchy kSZ power spectrum, onl y 
upper limits (jReichardt et all l2012t iZahn et all 12012ft . 
A detection of the patchy kSZ power spectrum will 
depend upon the ability to properly model contribu- 
tions from the thermal SZ power spectrum (which de- 
p en ds on uncertain intr a cluster medium astrophysics 
e.g. iBattadia et all l2010t iShaw et all 120101 : iTrac et al.l 
201lHBattaglia et al.ll2012al), the homogeneous kSZ (e.g . 



Ostriker fc Vishniad (l986t iJaffe fc Kamionkowskil H 998: 
Ma fc Fry! 120021 : iZhang et all 120041: IShaw et al.l l2012Tl. 
the thermal SZ - CIB cross spectrum (]Reichardt et all 
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^ Excluded by \ 
i Lya-K at z=6 N v 



_ _ Planck + luiC 2 

Planck + 0.5U.K 2 = 

_._ CMB-Pol + 0.5U.K 2 ~E 
X Fiducial Value = 




Fig. 8. — Projected constraints on z and A z for our fiducial model 
of EoR given the estimated error bars for Planck and CMBpol, and 
optimistic hypothetical error bars for ACT-pol and SPT-pol. The 
Ict (inner ellipses) and 2a (outer ellipses) constraints for Planck 
projected error bars on r with 1/iK 2 (blue dashed lines) and 
0.5/^K 2 (red lines) estimated error bars for ACT-pol and SPT-pol 
on D ^=3000 Patchy kSZ measurements. The black dot dashed lines 
are the lcr and 2cr constraints for CMBpol projected error bars on 
t and 0.5^tK 2 estimated error bars for ACT-pol and SPT-pol on 
D^_3ooo Patchy kSZ measurements. The grey region across the 
top is excluded by the zero transmi ssion of rest-fram e Lya flux at 
z > 6 in spectra of quasars (Lva-K: lFan et al.ll2006bh . 



|20H lAddison et ail [20121: iZahn et all [20121) . an d the 
infrared and radio sources (see iDunklev etaTI [2011 
iReichardt et al.ll2012L for these models). We choose two 
hypothetical error bar values for ACT-pol and SPT-pol 
measurements of -Df = 3ooo- 

The likelihood surface is C oc e~ x / 2 , and the lcr 
and 2cr contours are the 68% and 95% probability sur- 
face. Figure [8] show the forecasted constraints on z 
and A z from future measurements. Results from Planck 
and a detection of the patchy kSZ power at I = 3000 
will constrain z to ~ 5%. The value of A z begins to 
be constrained when we combine the upper limit from 
the opacity of the Lya forest (Lya-K) at z = 6 in 
quasar spectra. Here the upper limit from the opac- 
ity of the Lya forest is derived by converting our A z 
to z(x c = 50%) - z{x e = 90%). Figure |U illustrates that 
an experiment like CMBpol will tighten these constraints 
tremendously. While we focused our analysis on r and 
-0^=3000 it is possible that the low-^ EE power spectrum 
from CMBpol could provide constraints on A z as well, 
which would make our forecasted constraints even better. 

5. CONCLUSIONS 

Using a new semi-analytic model for constructing the 
reionization-rcdshift field from any density field, we 
made predictions for the \ow-£ polarization power spec- 
trum and the high-£ temperature spectrum measure- 
ments of the CMB. We demonstrated that combining 
measurements of the EE power spectrum with the patchy 
kSZ amplitude at £ = 3000 constrains both the mean 
reionization-redshift and the duration of reionization. 
The measured EE power spectrum from WMAP and the 
predicted improved spectrum from Planck will constrain 
z, but cannot discern between models for reionization 
with extreme durations. 



The shape and the amplitude of the patchy kSZ power 
spectrum depend on both the duration and the mean 
redshift of reionization. At I = 3000, where the ampli- 
tude of the kSZ power spectrum is currently constrained, 
we find the patchy kSZ power at £ = 3000 ranges from 
0.87 - 2.42 fiK 2 with a z = 10 (this z matches the current 
WMAP 7-year constraints). The largest kSZ signals cor- 
respond to long duration reionization models. We found 
a simple scaling law for the patchy kSZ power spectrum 
amplitude at £ = 3000 as a function of z and A z , which 
makes model fitting to observed spectra trivial. Using 
this scaling law and constraints from WMAP on the z 
and the lower limit from EDGES we place a lower limit 
on the patchy kSZ amplitude at £ = 3000 of ~ 0.4 fiK 2 . 

The amplitudes we find for the patchy kSZ power spec- 
tra are lower than previous model predictions. The 
differences between o ur work and the earl i er work on 
the patchy kSZ signal (jMcQuinn et al.ll2005t IZahn et all 
20051 llliev et~aT1 l2007) are that previous work used much 
smaller volumes to calculate this signal. When small 
volumes are used to calculate the patchy kSZ signal, we 
show that the amplitudes of the power spectra are biased 
to larger values. 

Many of our reionization models are con sistent with 
the ti ghtest upper limit constraint from SPT (IZahn et al.l 
I2012D . including our fiducial model, which has a z = 10 
and t = 0.085 to be consistent with the current WMAP 
7-year constraints. This consistency is achieved without 
the need to invoke more exotic or unphysical models for 
reionization that were necessary for the previous models 
of the patchy kSZ in order to fit within these constraints. 
In the event that the measured values of r and z decrease, 
the patchy kSZ power for all models would decrease fur- 
ther. 

It is clear that the future measurements from the 
Planck satellite of the EE power spectrum will tightly 
constrain the mean reionization-rcdshift and CMBpol 
has the potential to do even better. This leaves mea- 
surements of the patchy kSZ power spectrum from high 
resolution CMB observations to constrain the duration 
of the EoR. The current upper limits of the patchy 
kSZ amplitude at £ = 3000 from SPT range from 2.1- 
4.9 fiK 2 depending on the assumptions made about the 
correlation between the tSZ and CIB, but future de- 
tectio ns are projected to be ~ 1 /xK 2 ([Reichardt et al.l 
120121 ). In order to measure the patchy kSZ, we first need 
to understand £ (IReichardt et al.1 120121 : lAddison et al.l 
120121 : IZahn etaLl I2012I ). have a good understanding 
of the contributions from the homogenous kSZ (e.g . 
Ostriker fc Vishniad 119861: Uaffe k Kamionkowskj 119981: 
Ma k Frvl 120021 : IZhang et al.1 120041: IShaw et al.1 \20m . 



and how the astrophysical uncertainties of the homoge- 
nous kSZ models. A measurement of the patchy kSZ will 
tighten the constraints on z and A z greatly. 

Combining these CMB constraints with neutr al hy- 
drogen measurements (e.g lLa Plante et all I2012D . such 
as redshifted 21cm signal that originates from the hy- 
perfine transition of neutral hydrogen (e.) iScott k Reed 
119901 : iShaver et all I1999L iZaldarriaga et all 120041 ) . will 
constrain the mean redshift and duration of reionization 
further, or provide new issues for reionization models to 
tackle. 
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